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■ Abstract 
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C/2 . 

J>^, In studies on Marangoni type motion of particles the surface tension is often 

Qh| approximated as a linear function of temperature. For deformable particles in a 

linear external temperature gradient far from the reference point this approximation 

> '■ 

\^ . yields a negative surface tension which is physically unrealistic. It is shown that 

<N , 

H. Zhou and R. H. Davis (J. Colloid Interface Sci., 181, 60, (1996)) presented 

(N 

' calculation where the leading deformable drop moved into a region of negative 

surface tension. With respect numerical studies the restriction of the migration of 
O i two deformable drops is given in terms of the drift time. 

Ph. 

> ■ 

I The bulk fluid motion induced by an interface has been studied for over a century. 

I One of the most interesting phenomena is the capillary motion of particles through a 
viscous fluid. Young, Goldstein and Block (|l|) and later Bratukhin (j^) performed the 
first systematic study of the migration of bubbles and droplets. As noted in review (§) 
the capillary motion arises due to gradient of the surface tension 7 at the interface as a 
result of a non-uniform temperature or surfactant distribution in the surrounding media. 
The surface tension gradient results in a tangential stress on the interface which causes 



1 



the motion of the surrounding hquid by viscous traction. Then, the droplet or bubble 
will move in the direction of decreasing interfacial tension. It is necessary to note that 
the normal component of the capillary forces arising during the motion may deform the 
shape of a particle (|^). Young, Goldstein and Block (|l]) and others have shown that in 
the limit of high surface tension (undeformed spherical particle) its motion is controlled 
by surface tension gradients only. Note that the motion of a deformable particle also 
depends on the surface tension itself. 

If the particle moves with constant velocity the transformation of a laboratory co- 
ordinate system to a coordinate system moving with the particle frame will essentially 
simplify the solution @, (|), (||). Let us denote the particle coordinate system moving 
with the droplet velocity U by O' and the laboratory coordinate system by O respec- 
tively. We consider the coordinate transform from O to O' in the case of a drop moving 
in the uniform external temperature gradient Ae^ (^, see Fig.l. For an arbitrary point 
F we obtain, 

R = R' + Ut, Vi(R,t) = v;(R') + U, Ti(K,t) = T-(B:) + AUt, [1] 

where i = 1,2 correspond to the inner and outer liquid phase, respectively, V is the fluid 
velocity, T' denotes the difference between the temperature T in O and a undisturbed 
temperature AUt at the center of O', R is a radius vector which points from O to F and 
t is the time. 

In the limit of an infinitely large surface tension the normal stress boundary condition 
is not modified under the above transformation P). However, in the case of finite surface 
tension, this boundary condition requires special attention. Usually, 7 is assumed to be 
linearly dependent on temperature or on concentration is linearized (^), 



7(R,T) = 7o(ro)+ 



(Tm-To), [2] 

T=To 



where dj/dT is a constant and Tq and 70 correspond to the reference values of temper- 
ature and surface tension, respectively. Note that for many cases d'j/dT < 0. Due to 
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the transformation of T the surface tension 7(R, T) is also transformed in the moving 
coordinate system, 

7'(R', T') = 7o(ro) -^^Ut+^ (r'(R') - To) . [3] 

The surface tension 7' is time dependent now. Recall that the surface tension must be 
positive 

7 > 0, 7' > 0. [4] 

From 0] and [|] follows an upper bound of the drift distance Ut in system O or an upper 
bound of the time of particle migration in the moving system O'. 

Ignoring the above restrictions results in the appearance of a negative surface tension 
in the course of the particle migration in finite time and thus may lead to a physically 
unrealistic behavior of the particle. This restriction is relaxed in the case of the un- 
deformed drop (|T]), (|^) and (|^) where the normal stress boundary condition is always 
satisfied. However, this is not true when the surface tension has a finite value. We noted 
that in the literature on thermocapillary migration of drops and bubbles no attention 
was paid to this point. For example, Zhou and Davis (0) first considered the problem of 
axisymmetric thermocapillary migration of two deformable viscous drops . The authors 
assumed a linear dependence of surface tension on temperature. In terms of we have 

7(x,) = 70 + — (T(x,) - To(x^)) , [5] 

where T(xs) is the temperature at a point x^ on the interface and To(xr) is a reference 
temperature. In an attempt to obtain a solution which is independent of the choice of 
x^, Zhou and Davis fix x,. to be the intersection point of the axis of symmetry with the 
surface of the leading drop, see Fig.l and their Fig.l in (|^. It is important to note 
that this choice of x,. means a coordinate transform from the laboratory frame to the 
coordinate system moving with the leading droplet. Hence, the normal stress balance is 
modified. The other boundary conditions and the governing equations remain the same 
due to the linearity of Stokes and Laplace equations (|^). For more details see (H) and @. 



Zhou and Davis (J^) give for the dimensionless surface tension in the moving coordinate 
system 

7(x,) = l-gf(x,), [6] 

where 7 = 7/70 is the dimensionless surface tension, q = aA{—d'~f/dT) is the di- 
mensionless rate of change of the interfacial tension due to temperature variation, 
T(xs) = (T(xs) — To{'Xr))/{aA) is the dimensionless temperature difference and a is 
the radius of the first drop. It can readily be seen that Eq.|^ defines surface tension 
which is positive for any time or migration distances. As we showed before, the correct 
transformation of the linear approximation [^] leads to a negative surface tension in finite 
time. The previous conclusion that physically acceptable solutions must be restricted by 
migration time contradicts Eq. 0. 

Let us derive the correct form of the transformed surface tension in terms of The 
problem of the migration of two droplets is evolutionary and it must be accomplished 
by a kinematic condition applied on the droplets' surfaces. The transformation from the 
laboratory coordinate system to the particle coordinate system are given by (^: 

R = R' + /*' U(t) rft, V,(R,t) = V;(R') + U(t), [7] 

Jti 

Ti(R,t) = T-(B:) + A r\J{t)dt. [8] 



The migration velocity of the droplet now depends on time and therefore the migration 
distance Ut on the right hand side of is given as an integral term, 

7'(R',T') = 7o(To) - U(t) dt+^ (t'(R') - To) . [9] 

In terms of (|^ we have for the dimensionless surface tension 

7(x,) = 1 + ^ / U(t) dt-q f (x,). [10] 



The integral term in Eq. |T0[ changes the scenario of a numerical calculation. The surface 
tension changes with time and it is necessary to keep 7 positive. 
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We shall now proceed to estimate the time when the surface tension of some point 
Xj, on the leading drop will not satisfy For simplicity let us stay in the laboratory 
coordinate system, for 7 = we obtain the relation [| 

gf(x,) = l, [11] 

where T(xs) = (T(xs) — To(xo))/(aA) and xq is a reference point in system O. It is 
readily seen that the dimensionless length of a spatial frame is given by X = ^. Then 
the maximum transformation distance of the leading drop is the difference between X 
and the initial position x^. For the case of equal material parameters considered by (^) 
we have a = 1 and the surface separation distance on the axes is ~ 1. Hence the length 
of the drops' drift is also ~ 1. Let us assume that the lower bound of the velocities for 
moving deformable drops is the velocity of non-deformable drops. For slightly unequal 
drops and a large separation distance between their centers the velocities are nearly the 
same and equal to the Young-Bratukhin value of 0.133.. Following Eq. [12] in (|^ we 
normalize this value with 2/15 because for inner and outer liquids the viscosity and the 
thermal diffusivity are equal. From this normalization procedure we obtain that the 
migration velocities are ~ 1. As a result the critical value of the migration time is ~ 1. 

We developed a numerical code for solving the problem of the motion of two de- 
formable viscous drops in an external temperature gradient Restrictions [Q were 
considered in the laboratory coordinate system O. In Fig. 2 we plotted the evolution of 
the minimum separation distance d between the droplets' surfaces in time. We chose 
a = 1, a = 0.5, q = 0.2 in terms of (|^, where a is the droplets radii ratio. The dotted 
curve confines the physical region where [Q is satisfied. The curves 1 — 10 correspond to 
different initial separations. Curve 2 is in agreement with the results given by Zhou and 
Davis (see Fig. 4 in (|^) and with the asymptotics for the non-deformable drops (§). 

Our computations show that the patterns of drops deformations are similar to those 
described by (|^) but correspond to smaller separation distances d. Note that our analysis 

^On physical grounds this hmit corresponds to phase transition. 
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is restricted by ^ while the results of (|^ lie in the physically unrealistic region. For 
initially spherical drops and an initial separation distance d = 0.01 Fig. 3 depicts the 
series of drops' profiles corresponding to the points a, b and c in Fig. 2. 

The author wish to thank A. M. Leshansky and T. Loimer for helpful discussions. 
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Figures 





Figure 2: Evolution of the separation distances for two deformable viscous drops under 
a linear external temperature gradient as a function of the initial separation (j^). The 
values of the parameters are the same as in (0): a = 1, a = 0.5, g = 0.2. 
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Figure 3: Deformation patterns for the initial separation 0.01 and a = 1, a = 0.5, 
q — 0.2. Figures a, 6, c correspond to the respective points on curve 10 of the Fig. 2. 
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